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Within a dynamical description of nuclear fragmentation, based on the liquid-gas phase transition 
scenario, we explore the relation between neutron-proton density fluctuations and nuclear symmetry 
energy. We show that, along the fragmentation path, isovector fluctuations follow the evolution of 
the local density and approach an equilibrium value connected to the local symmetry energy. Higher 
density regions are characterized by smaller average asymmetry and narrower isotopic distributions. 
This dynamical analysis points out that fragment final state isospin fluctuations can probe the 
symmetry energy of the density domains from which fragments originate. 



PACS numbers: 25.70.Pq, 21.30.Fe, 24.60.-k, 05.10.Gg 

The dynamics and thermodynamics of complex sys- 
tems present general aspects, of interest in different do- 
mains of physics. A rather important issue is the iden- 
tification of the occurrence of phase transitions. This 
is relevant for many microscopic or mesoscopic systems, 
from metallic clusters to Bose condensates and nuclei [J- 
0]. In particular, the analysis of two-component systems 
has recently evidenced new interesting features [J-l7[. 

Under suitable conditions of density and temperature, 
the nuclear Equation of State (EoS) foresees the possi- 
bility of phase transitions from the liquid to the vapour 
phases, a scenario often evoked to explain the multifrag- 
mentation phenomenon (8l-[lo|. As a consequence of the 
two-component structure of nuclear matter, constitued 
by protons and neutrons, a crucial role is played by the 
low-density behavior of the isovector part of the interac- 
tion and the corresponding term in the nuclear EoS, the 
symmetry energy [11], on which many investigations are 



rgy [li] 

imi. 



concentrated [12Ml6j . We stress that this information is 
essential in the astrophysical context, for the understand- 
ing of the properties of compact objects such as neutron 
stars, which crust behaves as low-density asymmetric nu- 



clear matter [I7|,|l8|. Moreover, the density dependence 



of the symmetry energy affects the structure of exotic 
nuclei and th e ap pearance of new features involving the 
neutron skin [19|. 

A connection between the characteristics of clusters 
emerging from nuclear fragmentation and the symmetry 
energy has been proposed, in the framework of macro- 
scopic statistical models However it would be 



important to explore this issue within a full dynamical 
description of the fragmentation process. Here we under- 
take such a kind of study for systems facing low-density 
(spinodal) instabilities and first-order phase transitions 
We investigate the coupling between the develop- 
ment of neutron-proton density fluctuations (isovector 
fluctuations) , to which isotopic properties are connected, 
and the growth of unstable modes of the total density, 
leading to the formation of nuclear drops (fragments). 
Thus the aim of this work is to examine the behavior 
of isovector fluctuations in rapidly evolving systems, to 



probe their possible relation to the symmetry energy and 
its density dependence. 

Theoretically the evolution of complex systems can 
be described by a one-body transport equation with a 
fluctuating term, that incorporates the effects of the un- 
known many-body correlations^the so-called Boltzmann- 
Langevin equation (BLE) [2^ . We follow the approx- 
imate treatment to the BLE presented in Ref. |25| . the 
Stochastic Mean Field (SMF) model. We solve the fol- 
lowing equation for the time evolution of the semiclassical 
one-body distribution function f{r,p,t): 



df df du df - 



dUext ^ df_ 
dr dp' 



where U[p\ is the self-consistent mean-field potential, 
Icoiiif] is the average collision integral and Uext{i') rep- 
resents an external, stochastic field. The coordinates of 
isospin are not shown for brevity. Within such a frame- 
work, the effective nuclear potential U is derived from en- 
ergy functionals that usually contain a term proportional 
to P, the symmetry energy Esym,{pJ) / A = Csym{p)l'^ 
(with / = {pn — Pp)/ p and p, pn, Pp denoting total, neu- 
tron and proton densities, respectively). 

Let us consider the behavior of nuclear matter pre- 
pared with a uniform density distribution po and with a 
Fermi-Dirac momentum distribution corresponding to a 
specified temperature T. The system is confined within 
a cubic box, with side L = 19 fm, with periodic bound- 
ary conditions imposed. The linear response analysis al- 
lows one to get a first insight into the fiuctuation dy- 
namics. For two-component matter one can identify two 
types of independent modes of the phase-space density: 
isoscalar-like modes, where neutrons and protons oscil- 
late in phase, and isovector-like modes, with neutrons 
and protons oscillating out-of phase. In particular, in 
the case of symmetric nuclear matter, the two types of 
modes correspond to oscillations of Z'* = fn+fp (isoscalar 
modes) and of f" — fn — fp (isovector modes). Let us 
denote by f^{p,t) {q — s,v) the Fourier transform, with 
respect to r, of the difference Sf^ — ~ Jq, where /q is 
the system initial phase-space density. The equation of 
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FIG. 1: Panel (a): The quantity F^j,f, extracted from SMF 
simulations, for stable nuclear matter in several density con- 
ditions and at temperature T = 3 MeV (thick lines). Thin 
lines show the density dependence of the symmetry free en- 
ergy Fsym ■ Full lines: asy-stiff EoS. Dashed lines: asy-soft 
EoS. Panel (b): The quantity F^^s, evaluated at the freeze- 
out time, as a function of the local density for unstable sys- 
tems with initial density pi (full lines), p2 (dashed lines), ps 
(dotted lines). Asy-Stiff EoS. Panel (c): The same as in panel 
(b), for asy-soft EoS. In panels (b) and (c) thick gray (cyan) 
lines represent the same results shown in panel (a) as thick 
lines. 

motion for these Fourier coefficients follows readily from 

+ *k v}^ 'Qp,^ ^ Q, - ^-^kk V , (2) 

Here dU^ / dp'^ represents the appropriate Fourier compo- 
nent of the derivative of the effective field with respect to 
the density and J^l^{t) is the Fourier component of the 
external field. Furthermore, is the Fourier trans- 

form of the density fluctuation 5p''{r). Finally, since we 
will restrict our analysis to rather low temperatures, in 
Eq.(2) we have ignored the average collision term Icoii, 
since its effect is relatively small [26|. 

For stable modes, the equilibrium variance cr^ asso- 
ciated with the fluctuation p^ is linked to the physical 
quantities that characterize the response of the system 
to the action of the external force F^, see Eq.(2). Ac- 
cording to the fluctuation-dissipation theorem [27|, one 

can write: cr^ = T/F'i{k), where i^9(fc) = + 1/7V), 
with TV = Jd'P^- We notice that F« is noth- 

ing but the second derivative of the system free energy 
density with respect to the density p'^. Considering the 
inverse Fourier transform of p^ we obtain, for the equi- 
librium spatial density correlations, in a cell of volume 
AV: 

(jI^Hav) =< 6p'^{r)5p'i{r) >= 



(2^E-fc^k=^<l/^^(fc)>k, (3) 

where the average extends over all k modes. 

Focusing on isovector modes, the potential rep- 
resents the Fourier transform of the symmetry poten- 
tial UsyMr)] = CryMr')] ■ 9.[\r - r'|), 

where Cf°^ denotes the potential part of Cgym and the 
smearing function gc is introduced to account for the fi- 
nite range of the nuclear interaction. Thus we obtain: 
F^{k) = 2C!^M 9.ik)/po + ^ 2F,,,„(fc)/po. 

We note that the function Fsym{po) = Fsym{k — 0) 
simply coincides with the volume symmetry free energy, 
that at zero temperature reduces to the symmetry energy 
CsymiPo)- We can write: 

< 1/F^{k) >k= ^ < l/F,y^{k) >k^ (4) 

Hence we find that equilibrium fluctuations of the isovec- 
tor density can be connected to an "effective" symmetry 
free energy fj^j that, owing to the k dependence of the 
symmetry potential, is smaller than the free energy Fgym- 

In asymmetric matter, the findings discussed above 
still hold for isoscalar-like and isovector-like oscillations. 

Now let us go back to the full non- linear equations (1), 
that are solved numerically with the test particle method 
[28| . We have performed SMF calculations for nuclear 
matter prepared at initial temperature T = 3 MeV and 
in several density conditions. Here we also take account 
of fluctuations in the isovector channel, which were ne- 
glected in Refs.m, H^. Isovector fluctuations can be ex- 
tracted from the model by simply rescaling the variance 
by the number of test particles employed in the simula- 
tion [s^. 

We adopt momentum-independent effective interac- 
tions corresponding to a soft EoS, with compressibility 
modulus K = 200 AleV. The coefficient Csym gets a ki- 
netic contribution just from basic Pauli correlations and 
a potential part, Cf°^, from the isospin dependence of 
the interaction. For the local density (p) dependence 
of Cfyl^ we consider two representative parametriza- 
tions: one with a linearly increasing behaviour with den- 
sity (asy-stiff), CP°l^{p) = 90 p (MeV), and one with 
a kind of saturation above normal der isity (asy-soft), 
C^ymip) = P (238 - 1009 p) (MeV) [H [29[. We notice 
that at the temperature considered in the calculation, 
which is within the typical range observed in multifrag- 
mentation 10], the symmetry energy Csym is very close 
to Ysym- As smearing function we take a gaussian 
with width cr = 0.9 fm. With this choice, for nucler 
matter at saturation density (po = Psat = 0.145 fm^^), 
Eq.(4) gives F,^^^ = 0.7 F,j,™. 

Let us consider first, for the sake of simplicity, the 
case of symmetric matter (/ — 0). We first concen- 
trate on isovector fluctuations for uniform matter at 
rest, where equilibrium conditions are fulfilled. Thus, 
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FIG. 2: (Color online) The quantity 4A(Z/yl)^) (see text) is 
plotted as a function of the local density, for systems having 
initial density pi (full lines), p2 (dashed lines), pa (dot-dashed 
lines). Curves are shifted for a better visibility. Black, asy- 
stiff EoS; gray (red), asy-soft EoS. 



in order to avoid the development of volume instabil- 
ities at low density 



we switch-off in the calcula- 



tions the isoscalar part of the nuclear potential. Then 
we calculate the isovector fluctuation variance Upv —< 
{Spn{r) — Spp{r))'^ >, where the average is performed 
over cells of volume AV — 1 fm^. The effective sym- 
metry free energy can be extracted from the numerical 
variance exploiting Eqs.(3,4). This quantity is displayed 
in panel (a) of Fig.l as a function of the matter density, 
for the two parameterizations of the symmetry energy in- 
troduced above (thick lines, full for asy-stiff and dashed 
for asy-soft), and compared with the corresponding sym- 
metry free energy Fsym ■ The numerical results generally 
go with the analytical estimation discussed above: Owing 
to the k dependence of the symmetry potential, the ex- 
tracted F^ff is lower than the symmetry free energy, be- 
ing reduced by about 30% at saturation density, and ex- 
hibits a density dependence connected to the asy-stiffness 
of the effective interaction employed in the simulations. 

The evaluation of the equilibrium isovector fluctua- 
tions of stable matter can be used as a benchmark for the 
general and more interesting case where unstable systems 
are let evolve. Calculations have been performed taking, 
as initial density po, three values inside the spinodal re- 
gion: pi = 0.0245 fm^^, p2 = 2pi and p3 ~ 3pi. More- 
over, for each case, we have considered symmetric matter 
(system (1), /i = 0) and asymmetric matter (system (2), 
h = 0.142). 

Now the system may develop density fluctuations, so 
locally the density gets larger (density bumps, leading 
to fragments) or smaller (vapour) than the initial value 
[m . The separation between the two regimes is smooth, 
so that the local density p may vary between zero and 
values around the saturation density. Our analysis is 
performed at the "freeze-out" time t = 200 fm/c, when 
isoscalar density fluctuations saturate. At this time, the 
average density of the regions having p larger than po 
goes from 0.064 fm~^ (in the pi case) to 0.10 fm"^ (p2 
case) and 0.12 fm^^ (pa case). 



Our aim is to investigate the behavior of isovector fluc- 
tuations on the short time scale (the "freeze-out" time) 
associated with fragment formation. Isovector fluctua- 
tions are evaluated as a function of the local density 
inside the fragmenting system, looking at the variance 
of the isovector density p"" in cells having the same lo- 
cal density p. As a measure of the isovector variance 
CTp^, we consider the quantity F^^^ — (p T)/{2AVap^), 
that coincides with F^jj: if equilibrium is reached (see 
Eqs.(3,4)). Results for F^^^, obtained in the case of 
symmetric matter, are displayed in Fig.l as a function 
of the local density, for the three initial density values 
considered, see panels (b) and (c). Quite interestingly, 
isovector fluctuations follow the local value of the sym- 
metry energy independently of the initial conditions of 
the system. Indeed the three curves associated with the 
different initial densities (full, dashed and dotted lines for 
pi, p2 and p3, respectively) are rather close to each other 
and they are also close, for each given local density, to the 
equilibrium results discussed above (here plotted as thick 
gray (cyan) lines), thus locally i^^^^ « ^eff- These re- 
sults indicate that, as soon as density fluctuations start 
to develop, a quick rearrangement of isovector fluctua- 
tions takes place, so that the equilibrium value corre- 
sponding to the new actual local density is approached. 
Indeed isovector-like oscillations are characterized by a 
much shorter time scale, with respect to the growth of 
the unstable modes [sij . Thus important coupling effects 
between isoscalar and isovector oscillations are emerging 
from the solution of the full non-linear Eqs.(l). 

Calculations have also been performed for the asym- 
metric system (2), leading to results very close to the 
ones displayed in Fig.l. In the latter case one can also 
discuss the isospin distillation mechanism, that induces 
a deviation of the local asymmetry from the system ini- 
tial value [2§|. In particular, we consider the following 
density-dependent quantity, derived from the symmetric 
system (1) and the asymmetric system (2): A(Z/A)'^ = 
(zJA)1-(zJA)1, where (zJA), (with z = 1, 2) represents, 
for the system (i) the average proton fraction of cells hav- 
ing the same local density p. This quantity is displayed 
in Fig. 2 as a function of p. The different curves corre- 
spond to the two EoS (gray (red) lines for soft, black 
lines for stiff) and the three initial densities considered. 
As a general trend, we observe the well known behavior 
of asymmetric systems: The low-density regions become 
more neutron rich, while high density regions are more 
symmetric, just in connection with the density depen- 
dence of the symmetry energy coefficient Csym{p)- Here 
what is interesting to notice is that the distillation mech- 
anism goes together with the density-dependent behav- 
ior of the isovector variances described just above. As 
shown by Figs. 1-2, large density domains are associated 
with larger F'^^^ (i.e. smaller fluctuation width crp« /p) 
and smaller asymmetry, whereas low density regions are 
on average more asymmetric, but also more fluctuating. 
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FIG. 3: (Color online) The quantity ln{Y2/Yi) is plotted as 
a function of A'^, for the charges = 1 — 10, in the case of 
the systems with initial density pi. The stiff parametrization 
in considered. Lines are t o guide the eye. The inset shows 
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density. 



4:A{Z/Ay ■ F^ff, as a function of the local 



Let us move to study the probability Y{Z,N) to find, 
inside a volume V, a given number of protons and neu- 
trons, Z — pn,v V and N — ppy V. pn.y and ppy 
denote neutron and proton densities averaged over V ^ 
whose sum yields the density pv- Here we consider 
V = (5.5 fm)^. The quantity Y{Z,N) is proportional 
to the probability of getting, in the volume V, a specific 
variation of the isovector density p", with respect to the 
average p^: P(p") « exp - {p" - 'ff'f /{2ap.). Using 
the identity p^/pv = I — {N — Z)/A and considering 
the equilibrium amplitude of apv (see Eqs.(3,4)), one can 
write, for the yield ratio between systems (1) and (2): 

lniY,/Y,) « [{I-hf iI-l2f]A F:ff/T, (5) 

where F^ff and the average asymmetry li are functions 
of Pv, (in our case /i — 0). We notice that the ratio 
Y2/Y1 does not depend explicitly on the volume V . After 
some algebra, Eq.(5) can be rewritten as: 

IniY^lY,) « [(ll-ll){N+Z)-2{h-L2){N-Z)]F:ff/T, 

(6) 

Expressing li in terms of the average proton or neutron 
fraction, we finally get: ln{Y2/Yi) « aN + with: 



a{pv) = A^{ZlAYF:ff/T, 



P{pv) = A^iN/AYF^ff/T. 



(7) 



Thus we recover the standard isoscaling relations [20j, 
but with density-dependent coefficients a{pv) and /3(/9v), 
linked to the effective symmetry free energy ^ j . 

The behavior of the exponent a is illustrated in Fig. 3, 
where we plot the quantity ln{Y2/Yi) as a function of N , 



for the charges Z = 1 — 10. Inspite of the implicit density 
dependence of the isoscaling parameters, we note that the 
slope a is the same for all charges. This result follows 
from the opposite trend, shown by Figs. 1-2, of the two 
quantities A(Z/A)^ and F^^^, so that the product keeps 
almost constant (see the inset of Fig. 3). More precisely, 
the quantities /S.{Z/AY and A(iV/A)^ go approximately 
— Ip^po (sij . counterbalancing the density de- 



as 



dp 



pendence of F^^^. In the case of a linear behavior of F^^p 
i.e. close to the conditions of our stiff case, the isoscaling 
parameters, Eqs.(7), would be exactly constant. How- 
ever also in the soft case (not shown in the figure), the 
exponent a is roughly the same for all Z values (within 
7%.) Within our framework, the nearly constant value 
of a (or j3) inside the fragmenting system could be at 
the origin of the experimental observation of the same 
isoscaling parameter for the several products issued from 
nuclear reactions (2l| . which in principle may originate 
from different density regions and/or have different av- 
erage density. Then, knowing a (or /3) and the average 
asymmetry of a considered reaction product, Eqs.(7) give 
the corresponding effective symmetry energy of the den- 
sity region from which it emerges. In other words, this 
analysis allows one to probe the local symmetry energy 
of clusterized systems. It should be noticed that this 
provides a different information with respect to the ex- 
traction of the total symmetr y e r ierg y associated with 
clusterized low-density matter [IJ, |32| . 

To conclude, in this paper we have undertaken a dy- 
namical study of the disassembly of two-component un- 
stable systems, focusing on the coupling between the 
development of isoscalar and isovector density fluctu- 
ations. For nuclear systems, we have shown that the 
amplitude of isovector fluctuations follows the evolution 
of the local density and approaches, within time scales 
compatible with nuclear reactions at Fermi energies, the 
corresponding local equilibrium value, that is linked to 
the density-dependent symmetry free energy. Thus frag- 
ment isospin fluctuations and isoscaling parameters are 
related to the symmetry energy at the fragment forma- 
tion density. These results are relevant to experimen- 
tal isoscaling analyses aiming at extracting information 
on the symmetry energy, a topic of strongcurrent inter- 



est in nuclear physics and astrophysics [12|, |17|, |33|-|37|. 
Though secondary decay effects are expected to reduce 
the sensitivity of these observables to the specific shape 
of the symmetry energy [s^] , this analysis should still al- 
low one to probe the range of values spanned within the 
low-density conditions reached in nuclear fragmentation 
reactions. Finally, it should be noticed that our study 
is perfomed within the semi-classical approximation. It 
would be interesting to introduce quantum fluctuations 
and investigate their influence on the relation between 
isoscaling, isotopic distributions and symmetry energy. 
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